Gravitational Cherenkov Radiation from Extended Theories of Gravity 
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We linearize the field equations for higher order theories of gravity that contain scalar invariants 
other than the Ricci scalar. We find that besides a massless spin-2 field (the standard graviton), 
the theory contains also spin-0 and spin-2 massive modes with the latter being, in general, ghost 
modes. The rate at which such particles would emit gravitational Cherenkov radiation is calculated 
for some interesting physical cases. 
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Issues coming from Cosmology and Quantum Field 
Theory suggest to extend General Relativity (GR) in 
order to cure several shortcomings emerging from as- 
trophysical observations and fundamental physics. For 
example, problems in the early time cosmology led to 
the conclusion that the Standard Cosmological Model 
could be inadequate to describe the Universe at extreme 
regimes. In fact, GR does not work at the fundamental 
level, when one wants to achieve a full quantum descrip- 
tion of space-time (and then of gravity). On the other 
hand, Einstein gravity seems to present problems also at 
infrared scales since the so called dark sector (dark en- 
ergy and dark matter) could be cured either detecting 
new fundamental particles constituting the largest part 
of cosmic matter-energy or claiming for alternative or ex- 
tended approach to the gravitational interaction. Given 
these facts and the lack of a final self-consistent Quantum 
Gravity Theory, Extended Theories of Gravity (ETGs) 
have been pursued as part of a semi-classical scheme 
where GR and its positive results should be recovered [l| . 
The approach is based on corrections and enlargements 
of the Einstein theory adding higher-order curvature in- 
variants and minimally or non-minimally coupled scalar 
fields into dynamics which come out from effective ac- 
tions of unified theories formulated on curved space-time 

Beside fundamental physics motivations, these theo- 
ries have received a lot of interest in cosmology since 
they naturally exhibit inflationary behavior which can 
overcome the shortcomings of standard cosmology. The 
related cosmological models seem realistic and capable of 
coping with observations. ETGs play an interesting role 
also in describing the today observed universe. In fact, 
the good quality data of last decade has made it possible 
to shed new light on the cosmic effective picture 

From an astrophysical point of view, ETGs do not re- 
quire finding candidates for dark energy and dark matter 
at the fundamental level; the approach starts from taking 
into account only the observed ingredients (i.e. gravity, 
radiation and baryonic matter); it is in full agreement 
with the early spirit of a GR that could not act in the 
same way at all scales 04H|- 



At a fundamental level, the main features of ETGs is 
the emergence of new gravitational modes that can be 
roughly classified as tensor and scalar modes related to 
massive, massless and ghost particles In this view, 
GR is an exception including only massless tensor modes. 
The further modes could have interesting effects at ultra- 
violet and infra-red scales and could play an important 
role both in the Standard Model of Particles and in grav- 
itational radiation Q. In particular, they could be con- 
nected to the symmetry breaking in the high energy limit 
(and investigated at LHC @) and have a signature in 
the cosmological stochastic background of gravitational 
waves 0, 0, y| . 

In this paper we want to investigate possible 
Cherenkov-like radiation effects related to these 

further gravitational modes that could have observational 
effects. For our aim, assuming any curvature invariant 
other than the Ricci scalarH a generic action for the 
gravitational interaction is [l| 

A = J d 4 x^df(R,P,Q). (1) 

where 

P = R aP R al3 

Q = R aMS R at ^ 5 . (2) 

Varying with respect to the metric, one gets the field 
equations: 

FGfiv = \g^ {f-RF)- (g^D - V^ V )F 

— 2 (fpR"R a u + /Q Ra/3"iiiR a 2) 

-g^ a Vfi{fpR ap ) - a(fpR^) 
+2V Q V^ [f P R a {fi S b eta v) + 2f Q R^f) , 

(3) 



We restrict to fourth-order theories which give the main con- 
tributions in renormalization process Q, but the following con- 
siderations can be extended to any higher-order model involving 
generic powers of the D-operator and combinations of curvature 
invariants. 



2 



where we have set 

F _ df df df 

F =dR> fp =dP' fQ =dQ 1 () 

and □ = g al3 V Q V p is the d'Alembert operator0. 
Taking the trace of Eq. ([3]), we find: 

\ [2/ — RF - 2V a V b ((/ P + 2f Q )R a ?)+ 



-2(f P P + f Q Q)} 



(5) 



Expanding the third term on the r.h.s. of (|5|) and using 
the Bianchi identity G .a = 0, we get: 

n(F+^(f P + f Q )R\ = ix 

[2/ - i?F - 2R^V a V (f P + 2/ Q ) - i?D(/ P + 2/ Q ) 
-2(/pP + f Q Q)\ . (6) 

If we define 



and 



F + -(f P + f Q )R, 



dV 



(7) 



= r.h.s. 

dp 

of ([6]), we get a Klein-Gordon equation for the scalar field 



U(t> 



dV_ 
~dp~ 



(8) 



It is clear that the scalar field <fi, assumes the role of a 
field induced by the further degrees of freedom of ETGs. 
Obviously, (p is identically zero in GR. 

Let us linearize around the Minkowski background and 
then we assume 



(9) 
(10) 



Then, from Eq. ([7]), we get 



6<f> = 6F + 1(6 f P + Sf Q )R + |(/ P0 + f Q0 )SR , (11) 

where R = R{r)^) = and similarly / P0 = §p-\ n „„ 
which is either constant or zero. Note that the in- 
dex indicates quantities calculated with respect to the 



The notation T^^j = ^(T^^+T^^) denotes symmetrization with 
respect to the indices (fi, u). 



Minkowski metric, that means no deformation. 8R de- 
notes the first order perturbation of the Ricci scalar 
which, along with the perturbed parts of the Riemann 
and Ricci tensors, are given by: 



SRfjiypa — - (dpdvhfta + dadfihvp — d a d v hp P — d p dph va ) , 

SR^u = ^ {dadvh^ + d a dph a v - dpd v h - Uh^) , 
SR. = dpdyh^ - Uh , 

where h = rj^hp^. The first term of Eq. (fTTj) is 

dF dF dF 

SF =dR^ SR+ dP^ SP+ 0Q^ SQ - (12) 

However, since SP and SQ are second order, we get SF ~ 
Fpo SR and 



5$ = \Fro + -(/po + f Q0 )J 5R . (13) 

Finally, from Eq. ([6]) we get the Klein-Gordon equation 
for the scalar perturbation 8<p 

U5(j) = ]- ^ -8(f) - 

3 F,ro + 3 (/po + /qo) 

+ \5R af3 d a dp(f P0 + 2/ Q0 ) - ^6Rn(f P0 + 2/qq) 
= m 2 S(j) . 

(14) 

The last two terms in the first line are actually zero since 
the terms /po, /qo are constants and we have defined the 
scalar mass as m 2 . = — , .> f° — — - c — r. 
Perturbing the field equations ^ we get: 

Fo(6R„ v - ^uSR) = 

-( V p„n - d„d v )(5<f> - l(f P0 + f Q0 )SR) 

-Vp»d a dp(f P0 5R a(! ) - n(f P0 SR^) 
+2d a d p (f P0 5R a { / v) + 2/qo 5R\^) ■ 



(15) 

It is convenient to work in Fourier space where the follow- 
ing substitutions have to be operated: d^hpu — > ik 7 hp V 
and 'Oihpu — > — k 2 hp V . The above equation becomes 

Fq{8Rpv - -r]p U 5R) = 
2 

(Vpuk 2 - kpk u )(8(j) - -(/po + Iqq)8R) 
+ilp»k a k p (f Pa 8R a P) + k 2 (f Pa 6Rp„) 
-2k a k b (f P0 SR a ^8 b v) ) - 4k a kf,(fQo 5R a ^ v f) . 



(16) 
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We can rewrite the metric perturbation as 



h^v h^ u 2 ^Iv-v V^t^hf j 



(17) 



and impose the standard gauge conditions d^h^ u = 
and h = 0. The first of these conditions implies that 
k fJl h fJ,L ' = while the second gives 



(18) 



Inserting into the perturbed curvature quantities, we get 

5R^u = ^ (2k tl k v hf + k 2 r]^h f + k 2 h^ v ) , 
SR = 3k 2 h f , 



k a k p 5R a (iiu) p 



((k 4 r]^ - k 2 k^kv)hf + k 4 h^) 



(19) 



Substituting Eqs. (|17|) - (fT9|) into (fT6|) and after some al- 
gebra we get: 



(20) 



Defining hf = —j? we find the equation for the pertur- 
bations: 



k 2 



L spin2 



huv = 0, 



(21) 



The solutions of Eqs. (|2T|) and (|2"2l can be written in 
terms of plane waves 

h/w = e M „("p) ■ exp(ik a x a ) + c.c. (24) 
hf = e(j?) ■ exp(iq a x a ) + c.c. (25) 

where 



k a = (Um spin2 ,jf) 



j2 . . 4- ^2 
spvn 



(26) 



and where m sp i n 2 is zero (non-zero) in the case of mass- 
less (massive) spin-2 mode and the polarization tensors 
e M y("p) can be found in Ref. [l8| (see equations (21)- 
(23)). Eqs. (|2"Tj) and (|24|l . contain the equation and the 
solution for the standard gravitational waves of GR [l9| 
plus massive spin 2 termfj. Eqs. (|2"2"|) and (|2"5|) are re- 
spectively the equation and the solution for the massive 
scalar mode (see also @). 

The fact that the dispersion law for the modes of the 
massive field hf is not linear has to be emphasized. The 
velocity of every ordinary (arising from GR) mode h^ v 
is the light speed c, but the dispersion law (the second of 
Eq. ([2^|) ) for the modes of hf is that of a massive field 
which can be seen as a wave-packet. Also, the group- 
velocity of a wave-packet of hf , centered in ~0 , is 



vg 



(27) 



which is exactly the velocity of a massive particle with 
mass m and momentum p. From the second of Eqs. 
and Eq. (|2"T|) it is straightforward to obtain: 



where m 2 



is = - /P0+4/Q0 - From Ec L- O we get: 



Uhf=m 2 s h f . (22) 



From Eq. (|2"Tj) it is easy to see that we have a mod- 
ified dispersion relation which corresponds to a mass- 
less spin-2 field (k 2 = 0) and a massive spin-2 field 
k 2 = i , F ° , — = — m?„,v,o. To see better this point, 

2 JP0 + 2/Q0 Apirtz 

let us note that the propagator for h^ v can be rewritten 
as 



G(k) oc 



1 



fc 2 A; 2 + : 



(23) 



spin2 



Clearly the second term has the opposite sign, which in- 
dicates the presence of a ghost energy mode (see also 

EM3). 

As a "sanity check", we can see that for the Gauss- 
Bonnet term Cgb = Q — 4P + R 2 we have fpo = —4 and 
/go = 1. Then, Eq. ([2~Tj) simplifies to k^h^ = and, in 
this case, we have no negative energy modes as expected. 



vg 



This means that the speed of the wave-packet is 



V^ 1 - v g) uj - 



(28) 



(29) 



Summarizing these results, we can say that consider- 
ing ETGs (which we have generically assumed as ana- 
lytic functions of curvature invariants) more gravitational 
modes than the standard massless ones of GR have to be 
taken into account. In fact, we can note that there are 
two conditions for Eq. (|T4| that depend on the value of 
k 2 . In fact we have a k 2 = mode that corresponds to a 
massless spin-2 field with two independent polarizations 
plus a scalar mode, while if we have fc 2 / we have a 



Here +c.c. means plus the complex conjugate of the preceding 
term. 
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massive spin-2 ghost mode and there are five independent 
polarization tensors plus a scalar mode Q ■ In particular 
taking ~y} in the z direction, we have 6 polarizations d 
and the amplitude can be written as 



,(+) j 



=(x) 

''111/ 

- ef,„ 



\1V ' [If 

exp(ik a x a ) 



c.c 



(30) 



where the first two terms describe the standard polariza- 
tions of gravitational waves arising from GR, while the 
others are the massive fields arising from higher order 
theories. At this point we can calculate the gravitational 
radiation. To this end, we have to calculate the average 
energy momentum tensor of a plane waves, that, in terms 
of polarization tensors, can be expressed as 



16-kG 



e Xp *e Xp --\e x x \ 2 



(31) 



where we have denoted with V the complex conjugate 
and the brackets (■■) indicate the average process. Note 
that the formula (f3"Tj) have to be calculated for each po- 
larization. 

For our aim, it is convenient to write the polarization 
tensor e M „ in Eq. (|24[) explicitly in terms of the Fourier 
transform of energy momentum tensor T^ v of a system 
that emits gravitational radiation (20| : 



e MI /(x,w) 



4G 



2V(k,a;) - i?^T A A (k,cj) 



where 



T„ v (k,w) = J d 3 x% v {x', 



uj e 



(32) 



(33) 



Here we are supposing that the radiation is observed in 
the wave zone, that is at distances r ~ |x| much larger 
than the dimension R = |x'| moa! of the source. 
The conservation equation for T^ v (x, t) is 



d 
dxf 1 



T»„(x,t)=0. 



Applying this result to the following equation 



T^(x,t) = T flv (^.,uj)e~ 



+ c.c. , 



(34) 



(35) 



4 The polarizations are denned in our 3-space, not in a spacetime 
with extra dimensions. Each polarization mode is orthogonal 
to another and it is normalized e ( j„e fl " = 28. Note that these 
further modes are not traceless, in contrast to the ordinary plus 
and cross polarization modes of GR. 

5 Note that here, e M „ , is the sum of the polarization tensors arising 
from higher order gravity, as shown in Eq. OOt . 



gives 



O.r 



■T i „(x,w)-iwT°„(x,w) =0. 



(36) 



Multiplying with e lk ' x and integrating over x, we find 
that T^„(k, to) is subject to the algebraic relations 



fc |t T" >I/ (k > w) = > 



(37) 



where fc M is the vector k ~ lux, k° ~ uj. 

Now let us calculate the emitted power per unit solid 
angle along a direction x ~ - [2(| 



dE _ 

dn ~ 

From Eq. (|3"Tj) . we obtain 
dE _ r 2 (k-x)k° 



2-i ^ 



dQ IQixG 
and after some algebra, we get 
dE Geo 2 ~ 



(x,uj)e Xv (x,uj) - -\e x x \ 2 



(38) 



(39) 



T A "*(k, w)T Ai ,(k, uj) - - \T X A (k, uj)\ 2 



, (40) 

The problem is solved once we have calculated the Fourier 
transform (l33l) 



^ = 2G| dojuj 2 [T x »%k,u)T x „(k,u)-±\T x (k,uj)\ 2 }. 

(41) 

It is convenient to express this result in terms of the 
purely space like components of T Al/ (k, uj). From Eq. 
(|3"T)) . we have 

T i(k,w) = -WTjifcw) , 



T 00 (k, W ) = fc , fe^,(k, U ) 



(42) 



where k ~k uj ~ x. Using these results in Eq. (|4"Tj) gives 
dE _ Guj 2 

dfl 7T 

where the projection operator is (l9| 



— = ^-A ijtlm (k)T^(k,uj)T lm (k,uj) , (43) 



A-ij,lm\k) — SuSji IkjkfyiSn -\- ^kikjkik m -\- 

--ktflm + ^ijklk m + -Simkikj . (44) 

If the energy-momentum tensor is a sum of individual 
Fourier components as 

T^(x, t) = e l "*T^(x, u) + c.c (45) 



then the field /i„„ will look like a sum of plane waves. 
The gravitational energy-momentum tensor is then given 



5 



by a double sum over these Fourier components, but all 
cross-terms drop out when we average over a time interval 
which is long if compared with the longest "beat period", 
that is, the reciprocal of the shortest frequency difference. 
The energy is thus given by a sum of terms like (|4"3"|). 
one for each frequency in the source [2l[. Suppose, on 
the other hand, that the energy-momentum tensor is a 
Fourier integral. Then in the wave zone looks like 
an integral over ui of the individual plane waves, and 
the gravitational energy-momentum tensor is given by 
a double integral J J dujduj' of products of these terms. 
The integrand, again, has time dependence exp{—i(ui — 
w')t), but now there is no longest beat period, so instead 
of computing the average power, we calculate the total 
emitted energy. This is given by integrating the energy 
over the whole time, and the effects is to replace the 
factors e~ tult e ZUJ 1 in the double integral for the power 
with 



exp(-i(uj - U)')t)dt = 2tt5(uj - w') . (46) 



The energy per solid angle, emitted in a direction k is 
thus a single integral: 



dE 

In 



2G / ui 1 



A -(k, W )T Al ,(k, W )--|T A A (k, W )| 2 



dui , 
(47) 



For a free particle, we have 



= m / dr ^ t S 4 (x - £(r)) (48) 



(49) 

where 7 is the Lorentz factor. This gives us the Fourier 
transform, 

Tfj,u (k, u) = m <J(<yw(l - k • v)). 

The argument of this 5-function is the Lorentz invariant 
quantity £ M . It tells us that the radiation goes out 
along a cone, just as with familiar Cherenkov radiation; 
and this can only happen for velocities v that are greater 
than c = 1, in our units. Now we put ([50)) into (|4"Tj) and 
get, 

^ = 2G fduju: 2 ^J( 7W (l _ k . v )) 5(0). (50) 

It is worth noticing that 

S (p) = J- / dre^> = ^L = £L. (51) 

where At is the time interval over which we observe this 
radiation process [T3| ■ 

Now we integrate over all angles and get the rate of 
energy emission, that is 



— ; — = G — / du> UJ. 

At 7 2 v 



(52) 



We need to introduce some cut-off for the integral over 
lu, and to this end, let us consider the basic quantum 
relation, E = huj, where the emitted quantum energy 
cannot be larger than the total energy E of the particle. 
Thus we get 



AE G rrf_ /£ N 2 



(53) 



Noting that E = 7717, and re-introducing the constant c, 
we get 



AE 
~At 



G m A c i 
2 h 2 v ' 



(54) 



By substituting the values of the constants, we get the 
order of magnitude of the process, that is 



AE 
~At 



2\ 4 

mc \ ( c 



eV 



-)x 10- 41 — . (55) 
v 1 sec 



At this point, some consideration are in order for 
Eq. ([54|) . We find that the emitted energy depends on 
massless spin-2 (standard graviton), massive spin-0 and 
massive spin-2 modes. Then according to the following 
values of the masses 



7^ massive 
mspin2 = massless 
m s 7^ scalar 



(56) 



we have different values for the energy. As an example, 
we can consider, the upper bounds on the graviton mass, 
m g , that come from direct or indirect observations of 
gravitational waves. A constraint on the graviton mass 
comes from indirect evidence for the emission of gravita- 
tional waves from binary pulsars [22| . The upper limit is 
of the order 



m Q - 7.6 x lQ- 20 eV, 



(57) 



Inserting this value in Eq. ([55]) . we obtain the emitted 
energy of the order 



M „ 3.33 x 10— 1 
At 



(58) 



Furthermore, graviton masses of the order m g = 10~ eV 
could be detected by observing the characteristic signa- 
ture of a a strong monochromatic signal in gravitational 
wave detectors due to relic gravitons at a frequency which 
falls in the range for both of space based (LISA) and earth 
based (LIGO- VIRGO) gravitational antennas (23[, in the 
frequency interval lO~ A Hz < f < 10kHz. In this case, 
we get 



On the other hand, considering a gravitational state of 
mass order m g = 100 GeV, that is a Higgs-like particle, 
we obtain 

and finally if m g = 1 TeV, we have 



In the last two cases, we are considering that massive 
gravitational aggregates could emerge in high-energy ex- 
periments like those now running at LHC (CERN). See 
also [|[ for details. 

These results are a clear indication that a gravitational 
Cherenkov radiation could be detected ranging from very 
low energy scales up to TeV scales as soon as further grav- 
itational degrees of freedom are considered. A detailed 
experimental analysis will be reported in a forthcoming 
paper. 
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